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Abstract 

Recently the radiative B decay to the strange axial- vector mesons, B — )■ J^'i (1270)7, has 
been observed with rather large branching ratio. This process is particularly interesting 
as the subsequent Ki decay into its three body final state allows us to determine the 
polarization of the photon, which is mostly left- (right-) handed for B{B) in the SM 
while various new physics models predict additional right- (left-)handed components. A 
new method is proposed to determine the polarization, exploiting the full Dalitz plot 
distribution, which seems to reduce significantly the statistical errors. This polarization 
measurement requires however a detailed knowledge of the Ki — Knir strong interaction 
decays, namely, the various partial wave amplitudes into the several possible quasi two- 
body channels, as well as their relative phases. The pattern of partial waves is especially 
complex for the i^i(1270). We attempt to obtain the information through the combination 
of an experimental input and a theoretical one, provided by the ^Pq quark-pair-creation 
model. 



1 Introduction 



The 6 — 57 process has been playing important roles to understand the electro-weak in- 
teraction of the Standard Model (SM). The Glashow-Iliopoulos-Maiani mechanism shows 
that in the SM, the fiavour changing neutral current such as 6 — S7 is forbidden at the 
tree level but only occurs through a loop level diagram. Inside of the loop, heavy parti- 
cles, much heavier than b quark, can propagate. Therefore, the 6 — > S7 process can be 



1 



used to probe indirectly such heavy particles, namely top quarks in the case of SM or yet 
unknown particles introduced by given models beyond the SM. 

By now, the branching ratio of the inclusive B — > Xg'-f process is measured at quite a 
high precision {Br{B X,7)exp = (3.55 ± 0.24 ± 0.09) x 10"^ p. The SM theoretical 
predictions for this observable is obtained at the next-to-next-to-leading order in QCD 
[Br{B — )■ Xs'j)th = (3.15 ± 0.23) x 10~^ [2]) and they are relatively in good agreement 
with the experimental value. However, these predictions have theoretical uncertainties 
coming from the CKM matrix element as well as various kinds of QCD corrections. As a 
result, even if we add some new physics contributions to the theoretical predictions, the 
total branching ratio often agrees with the experimental value within those theoretical 
uncertainties. While tremendous efforts in order to improve the precision of the theoretical 
prediction have been made so as to match to the experimental precision, which could 
become even higher in the future machines, it is necessary to investigate the characteristics 
of the particles inside of the loop of the 6 — )■ 57 process using another kind of observables. 
In this article, we discuss a measurement of the circular-polarization of the photon of the 
6 — 57 process, which the left- and right-handedness of the couplings of the interactions 
among the particles inside of the loop. In the SM, the fact that the W boson couples 
predominantly to the left-handed quarks induce the photon polarization to be mostly 
left-handed. On the other hand, many new physics models contain new particles which 
couple differently from the SM. Therefore, the measurement of the photon polarization 
can be an useful tool to distinguish the interactions of the particles inside of the b sj 
loop from the SM-like one. 

Although, there have been several proposals for how to measure this photon polar- 
ization, its precise measurement has not been achieved yet. In this paper, we revisit 
the method proposed by Gronau et al. [3] (the GGPR method in the following) us- 
ing the exclusive B — )■ A'rcsonancc7 followed by the three body decay of the A'resonance- 
Most interestingly, the Belle collaboration recently observed one of these decay channels, 
B A'l (1270)7 ~^ (A'-7r7r)7, and found a relatively large branching ratio: Br{B~^ — )■ 
K^{1270)-f) = (4.3 ± 0.9(stat) ± 0.9(syst)) x 10^^ g], which dominates over the decay to 
A'l (1400), previously studied in detail by GGPR [3]. Thus, it is interesting to reconsider 
the feasibility of this method. In this article, we introduce a new variable, u, which is 
originally proposed by Davier et al [5] for the r polarization measurement at LEP (the 
DDLR method in the following). As we show later-on, the fact that the decay width of 
B — A'i7 — )■ {Ktth)'-/ process depends only linearly on the polarization parameter 
allows us to use the variable u in our study. And the simplification of the fit by using 
u makes it easier to include to the fit not only the angular dependence of the polariza- 
tion parameter but also the three body Dalitz variable dependence, which improves the 
sensitivity to the polarization parameter as also pointed out in [5]. On the other hand, 
the new radiative decay, to A'l (1270)7 instead of the A'l (1400)7, implies a more complex 
pattern of hadronic decay channels, not only through K*7r, but also through Kp and a 
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possible KTi. In this work, we discuss, in details, the hadronic parameters required in this 
analysis. In particular, having various difficulties to extract them fully from the currently 
available experimental data, we attempt to evaluate them with a help of the so-called ^Pq 
decay model. 

In section 2, we show a demonstration of the photon of the b s'-f being predominantly 
left-handed in the SM. We also discuss briefly the contamination from the right-handed 
polarization. In section 3, we derive the master formula for the decay width of the 
B — )■ Ki'-f — )■ {K'ir7r)'y and the hadronic parameters needed in this formula are evaluated 
in section 4. In section 5, we introduce the new variable u. We show our numerical results 
in section 6, including the comparison of the sensitivity of the DDLR method with the 
other proposed methods . Section 7 is our conclusions. 



2 Photon polarization of the 6 — > 57 in SM 

In SM, the quark level 6 — ?■ 57 vertex without any QCD correction is given as: 

where q = Pb — Ps with pi, and Ps four-momentum of b and s quark, respectively, F2 is 
the loop function, whose expression can found in [6]. When we fix the three momentum 
direction, namely the q direction as +z in the b quark rest frame, one can compute 
explicitly the helicity amplitude and we readily find that the ffist (second) term is non- 
zero only when we multiply the left (right)-handed circular-polarization vector, which is 
defined as: ^ 

e^ = -^(0,l,-z,0), 4 = — (0,l,z,0) (2) 

Since ms/rrib — 0.02 <^ 1, the photon in 6 —> 57 in SM is known to be predominantly 
left-handed. 

Once we include the QCD corrections, the other types of Dirac structure contribute 
and the above conclusion can be slightly modified. The result can typically described in 
term of the following effective Hamiltonian: 

Heff = -^v*Vtt I V + cmOjM + cmom] (3) 



V2 



J=l 



where Cj are the short-distance Wilson coefficients that can be calculated in perturbation 
theory and Oi are the local four-quark operators {i = 1 ... 6) and and Ogg are 
the electro-magnetic and chromo-magnetic penguin operators, respectively. The /i is the 
renormalization scale which is chosen as around m^. Note that Oj^ = j^mi,SLa(^^i,F^''bRa 
is equivalent to the first term in Eq. In addition to the small rris/mi, contribution. 



*The term proportional to (the second term in Eq. ([TJ) is neglected in this expression due to its 
smallness. 
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there is potentially non-negligible right-handed pollution due to the perturbative and 
non-perturbative contributions. A numerical estimate for them is extremely important 
while it is currently not available for B — t- Ki'-f. On the other hand, many efforts have 
been made in the case of i? — )■ K*'y using various Q CD-based approaches [7]-[TT]. Note 
that the time-dependent CP asymmetry of -B — ?■ K*'-/ can also be used to determine the 
photon polarization (see Section 6 for some discussions). To have an idea, it is found that 
the most recent estimate for B — )■ K*'~f [TT] shows that the right-handed correction is less 
than 1%, while another estimate [U] shows that it can be up to 10%. 

On the other hand, when we consider the new physics contributions, the right-handed 
contribution can be significantly enlarged from different types of Dirac structure that 
those new physics models can induce. It should be emphasized that there are many 
of the new physics models which can accommodate e.g. a large coefficient to the right- 
handed electro-magnetic operator (Oy^ with the subscripts L and R interchanged) without 
contradicting to the precise measurement of the inclusive B — > Xg'j branching ratio, for 
example Everett et al. [12]. 

3 The B Ki^ [Ktttt)^ decays 

3.1 Master formula for B — )► Ki^ — (PiP2^3)7 decays 

Due to the angular momentum conservation and the fact that 5-meson is a pseudo- 
scalar meson, helicity is conserved. Thus in order to determine the photon polarization 
it is sufficient to measure the polarization of the axial- vector meson Ki{l^) through its 
three body decay. As the physical final state Ki must have either left- or right-handed 
polarization, the decay width can be written as 

r(B ^ Zi7) = r(B ^ k,lil) + r(5 ^ ¥1^7^) (4) 

If we assume the narrow width of Ki, one can write the total quasi- four body decay width 
by these two terms, respectively, followed by the three body decay widths 

r(K,L ^ P1P2P3), r(K,R ^ P1P2P3) 

However, the width of the i^i's is not really negligible (r(ii'i(1270)) = 90 MeV, 
r(Jri(1400)) = 174 MeV according to PDG). Therefore, we present for completeness, in 
the following, a prescription that includes the initial state width of the Ki decay into 
the three body final state assuming the Breit-Wigner form, but which will not be used 
in practice. The Breit-Wigner factor is common to both polarizations and appears in 
modulus squared (therefore, its phase does not affect the crucial interference between the 
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J and J* terms below). Thus, our decay widths can be written as: 

rfr(5^Ki7-^(PiP2P3)7) 



dsdsi3ds23d COS 9 



oc ^ T{B ^ -ft'ipoi.Tpoi. 

pol. =L,R 



X 



dr(i^ipoi, ^ P1P2P3) 



(5) 



X 



dsdsi2,ds2'id cos 6' {s — m)^^ ^ 

where s = {pi + P2+ PsY is the off-shell "p^" of the Ki and Sij = {pi + pj)^ with pi to 
be the four-momentum of the final state Pj. Defining the —z direction as the photon 
direction in the Ki rest frame (see Fig. [T]), the 9 is given as cos 9 = ' P'^^P'^ 



IP1XP2I 



7^ {P2) i 



Pl XP2 




7r+(j9i) 



Figure 1: The Ki — )■ Kirir decay plane in the rest frame of Ki. Defining the —z direction 
as the photon direction, the 9 is given as cos ^ = f |^|^f?2| 



The kinematic distribution of this three body decays carries the information of the 
Ki polarization. It is important to notice that the polarization information we would 
like to obtain is the difference between r(P — Kil'Jl) and r(P — )■ Ki^'Jr) in Eq. (|5]) 
while experimentally, only the L.H.S of this formula, i.e. the total decay width can be 
measured. Thus, the high sensitivity to the polarization information can be achieved only 
if there is a significant difference in the decay distributions between Kn and Km. 

The differential decay width of Kh^r decay can be described by the helicity amplitude, 
J'fj_, which we define as: 

M(K,L,R ^ P1P2P3) = 4,L,«^M- (6) 

Considering that J"^ represents the decay amplitude of the Ki decaying into three pseu- 
doscalar mesons, we can parameterize it in terms of two functions Ci^2'- 

Jfj, = Si3, S23)piti - ^2(5, Si3, S23)p2ii (7) 
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where we omitted to write explicitly the Dalitz and angular variable dependences of J^. 
Note the s dependence of the coefficients, which means that in principle there could be 
some dependence on the off-shell of the Ki. Nevertheless, this dependence is not 
important as soon as the integration is limited to the Ki bump, especially for the ratio uj 
which is the relevant quantity in our method (see next section). The detailed expressions 
of Ci2(s 

; ■Sis, ■523) for given channels are derived in the next section but here we note 
that Ci^2('5, Si3, S23) can contain complex numbers. Using the definition of the helicity in 
Eq. ([2]), one can easily find in the Ki reference frame 



dsdsizds2zd COS 6 



CX |7W(J^iL,R->PlP2P3)r 

(8) 

oc ^|jfp(l + cos'^) T hm[n- {J X J*)] cose 



where n = F}^^-3. so that: 

blXP2| 



= |Cini;i|2 + |C2npir-(CiC2* + c*C2)(pi-P2) (9) 

n-{jxj*) = -{CiC; - ClC2)\pi X p2\ (10) 

where pi ■ pj = EiEj — {sij — mf — m^)/2 and \pi x pj\ = Pi ■ Pjtan~^0 with Ei = 
{s - Sjk + m^)/{sy^) and = cos^^[{pi ■ pj)/{\pi\\pj)]. 

It is worth mentioning that the difference between the left- and right-handed polar- 
ization amplitudes comes from the second term of Eq. which, to be non vanishing, 
requires the amplitude to contain more than one amplitude with non-vanishing rela- 
tive phase. Such a condition can be nicely realised in this decay channel since when Ki 
decays into three body final states through more than one intermediate two-body chan- 
nels, such as K*7T and Kp, there is the non-vanishing relative phase originated from their 
Breit-Wigner forms (based on the isobar model). 

Finally, the master formula is obtained in terms of the polarization parameter (the 
B ^ Ki form factor is skipped as being a common factor) 

rfr(5^Zi7->(PiP2P3)7) 1 

oc 



with 



dsdsi^ds2zd cos 6 {s - m\^Y + m\V\^ ^^^^ 

X ^\J? {I + cos^ 9) + \^hm[n ■ {J X J*)] cos 

^ r(p^i^i7) ^ ' 

which agrees with the expression in [3]. 

Since the two Ki resonances, Ki{127{]) and -ft'i(1400), are rather close each other, one 
could expect interference between them, which should then be taken into account in the 
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formulas. On the other hand, the Belle data |1] shows a suppression of the i^i(1400). 
In [131 E] it has been shown that such a suppression can be explained by taking into 
account the fact that these two states are the mixture of l^Pi and l^Pi sates and a 
reasonable choice of the mixing angle can explain such a suppression (see section 4 for 
more detailed discussion on the mixing angle). 

3.2 The Ci 2 functions for the Ki — > Rtttt decays 

In this section, we derive the Ci^2 function, which is defined in Eqs. and ([7]) for the 
1270/ 1400) decay. The three body decay channels of the i^i (1270/1400) are the ttttK 
final states. We first assume that this three pseudoscalar meson final state comes form 
the quasi- two-body decay through a vector meson, namely p or K*. The different decay 
channels and the possible vector resonances for ii']^(1270/1400) and i^'°(1270/1400) are 
listed below. 



I : 


ir+ (1270/1400) 






(13) 






P+ 






II : 


i^i+(1270/1400) 




)7T+{p,)K + {p,) 


(14) 






















III : 


i^?(1270/1400) 




7r-{p2)K + ip,) 


(15) 












IV : 


/s:°(1270/1400) 


^ 7r+(j9i 


K-+ 


(16) 



K*+ 

The decay amplitudes for these decay channels can be written as the sum of the amplitude 
with different intermediate vector meson channel: 

V 

where Pi, 2,3 represent the final state mesons carrying the momentum ^1,2,3 as assigned in 
the Eqs. (IT^ to (ITB]) and V represents the vector meson resonance. The Clebsch-Gordan 
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coefficients, Cijk, for each intermediate channel are given as: 

M,v{Kt ^ir+{p,)ix-{p,)K\p,)) = -^Mf;;p^)p^ - ^M^l^p^^p^ 

Using the detailed expression for the quasi- two-body decay amplitude Aijp^p^^p^ given 
in Appendix A, we find: 



M{K^L,R ^ 7r7rK)A=i^jv = e^^^^pj^ (18) 

= (S, Si3, S23)PlM - C^is, Si3, S23)P2/. (19) 



with 



(20) 



where 



^5 = 9vp.p,'BWv{s^,)[-fv + hvV^{E, - Ej) - A, 



(21) 



with A,, = + ^v^(^« + ^J-)]- 



4 Hadronic parameters and their estimation in the 
^Pq model 

The next step is to obtain the coupling constants and the form factors determining the 
above functions Ci^2, i-e. the following hadronic parameters 

Qpn-TT, gK'-KTT, fv, hv 

Noting that there are total of four fy and hy (V = p, K*) for each i^i(1270) and i\:i(1400), 
we have ten free parameters in this decay mode. One may consider the relative phases 
between the form factors fv and hy, which increases the number of free parameter. 
However, this phase could actually be determined theoretically or experimentally. 
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Ideally, these parameters should be extracted from the same experimental data of the 
B —J- i^i7 decay, which allows us a hadronization model independent analysis. For the 
present moment, one can use the other experimental information such as p — )■ tttt and 
K* — )■ Kir for the g couplings and Ki — )■ Kp and Ki — K*'it for fv and hy form factors. 
We first present in the following subsection how to relate those experimental information 
to our hadronic parameters. 

4.1 Hadronic parameters 

The VPP coupling constant: gvpp'- 

The gvpp coupling can be extracted from the partial decay width of the vector mesons. 
These are well measured for V = K* so that we can obtain this coupling rather precisely. 
The partial decay width can be written as: 

nv ^ P,P,) . (22) 

where |p| = ^y {my — (mi + m2)^)(m^ — (mi — m2Y)/2my. Then, using the experimen- 
tal values of p and K* widths, we find 

gp.. = (5.98 ± 0.02), gK^Kn = (5.68 ± 0.05) (23) 



The Ki — > VP form factors fy and hy 

To describe the Ki — )■ VP decay, we used two independent form factors fy and hy 

{V{py,e^^^)PM)\^HA\K,{pK,,e^-^^)) = ef^^T^^e^P* , T^'^ = fyg^'' + hyp^yp^^^. 

(24) 

On the other hand, the Ki — )■ VP can also be written in terms of the helicity amplitudes 
for the two possible +z spin projection of Ki and the vector meson, {Xki,^v) = (0,0) 
and (1,1). These two amplitudes actually can be written in terms of common partial 
wave amplitudes. Thus, when we expand them up to L = 2, we can equivalently write 
these helicity amplitudes by the two partial wave amplitudes |15j : 

{V{-pk,Xy)Pk{Pk)\AHA\K,{6,XKj) = {A^y + V5{2,0;l,Xy\l,Xy)A^)Dll^^^^{Qy) (25) 

where Ay^ are the partial wave amplitudes. Then, these amplitudes can be experimen- 
tally extracted through the partial wave analysis of the Ki — >■ VP processes using: 



r{Ki^VP)s-...o = ^\Ay 

ousa 



T{K,^VP)d-...c = P^K\' 

OUSa 
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\PV\ I AS 12 

(26) 



Comparing to Eq. (124|) and Eq. (l25l) . we can immediately find the relation between the 
two form factors and the partial wave amplitudes {fy, hy depend in general on Ey: 



72^ 



V 



Ev 



4^ 



1 + 



2El 



./iy 



(27) 



Partial wave analysis of Ki VP process has indeed been performed by the ACC- 
MOR collaboration [16] and a very precious information related to Ki meson has been 
extracted, which constitute the basis of the PDG entries. It is the currently available 
most extensive study of the Kmr channels, with full angular distributions analysis, deter- 
mination of relative phases between all amplitudes. On the other hand, the interpretation 
of the ACCMOR data contains various problems in the theoretical point of view, or even 
empirically 0. We will come back to some of these issues later in this section. In any case, 
we found that it is currently impossible to extract all the parameters from experimental 
data. Thus, we need a help of theoretical model inputs for this reason. In the following, 
we try to use the so-called ^Pq model, that is an intuitive model describing the decay by 
the creation of a quark- ant iquark pair. 

^We give some examples of those problems: 

• The resonance study is done by using the K-matrix method. Therefore, to match the information 
obtained from their analysis to our Brcit-Wigner parametrisation is not a simple task. Masses and 
widths must be recalculated. In addition, the authors use a complex phase space. 

• ACCMOR results are obtained by using particular models for strong interaction production 
through Kp scattering, like the Deck effect. Moreover, and more worrying, the D wave in 
i<ri(1270) — >■ K*Tr is depending strongly on the production transfer t. This fact may escape 
attention of PDG readers, because it averages between the two sets of data (high t, low t). 

• There arc only limited amount of data on the D-waves. They are poorly measured in K*'!t. For 
the D-wave amplitude i4ri(1270) — >■ Kp, there is no information... 

. Two other items are found to be important issues for the determination of polarization, the question of 
relative phases between the various partial waves, and the ktt channel; we devote to them two separate 
paragraphs below. It is to be noted that the Babar collaboration [T^ has performed a reanalysis of the 
ACCMOR data ; it contains useful complementary information, with somewhat different results for the 
parameters. On the other hand. rcf. [18] is a new, completely independent, analysis, which comes to 
certain conclusions differing from ACCMOR, especially for the ktt channel. 
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4.2 Estimating the hadronic parameters in the ^Pq model. 

The ^Pq model @ has the advantage that it provides rather complete predictions-in par- 
ticular the model fixes the coupling signs, ratios of K*7r to Kp couplings, D/S ratios 
and the full set of couplings once the quark pair-creation constant 7 is fixed. The model 
fixes the ratio of two independent couplings (C = ±1) which is left free by the SU{3) 
symmetry. Another illustrative example of its specific usefulness is the prediction of a 
very small decay of ifi(1270) into -R'g(1430). On the other hand, we must stress that it 
is a very approximate model, not claiming to be always quantitative. Its main drawback 
is that it is essentially non relativistic. 

There are two independent Ki states in the quark model, the l^Pi and l^Pi states, 
which are called Kia and Kib respectively. With the ^Pq model, we can predict the 
decay rates of these two states. However, these are not the physical mass eigenstates 
A'i(1270, 1400): it has been known that the observed hierarchy of decays into K*'ir and 
Kp, can be nicely explained by considering that the physical states are a mixture of Kia 
and KiB with a mixing angle 9ki '■ 

|iri(1270)) = \Kia) sinOK, + \K,b) cosOk, 

(z8) 

\K,{UOO)) = \K,a) cos 9k,- \K,B)sm9K, 

Then, the 8 independent amplitude A^.^^) each for A'i(1270) and i^'i(1400) can 

be reduced to the 4 amplitudes A^^, and one mixing angle 6k,: 



^K,il270)^K*n/Kp = S K* / p{V2 Sm 9 K, T COS 9k,) 

^K,{i270)^K'^/Kp = Dk*/p{- sin 9k, tV2 cos 9k,) ^^^^ 
'^K,{iAoo)-^K*7T/Kp = SK*/piV2 COS 9k, ± siu 6'a'i ) 

AK,{1400)^K*n/Kp = D K* / p{- COS 9 K, ±V2sm9K,) 

The Sk*/pi Dk*/p amplitudes are expressed in terms of the hadron wave functions and of 
the quark-pair creation constant 7 in this model. Having these model parameters fixed, 
we can obtain the mixing angle by fitting to the experimental data of Ki decays. The 



■i-This model has been first developed in [TO] and then extensively discussed by the group around 
N. Isgur [20l[21][22] in Canada. It has been already used by Blundcll et al. in the present context [23] . 



11 



available experimental information is listed below: 



5r(i^i(1270) ^ K*tt) : 5r(A'i(1270) ^ Kp) = (16 ± 5)% : (42 ± 6)% 
Br{Ki{UOO) K*Tx) : Br{Ki{Um) Kp) = (94 ± 6)% : (3.0 ± 3.0)% 



Sr(i^i(1400) ^K*7r)5-wavc 
Er(i^i(1270) ^ A-*7r)^_^ave 
5r(Ki(1270) ^ K*7r)5-wave 

Br{Ki{mO) ^ K*vr)5-wavc 



0.04 ±0.01 
0.54 ±0.15 (low t) 
2.67 ±0.95 (high t) 



(30) 



where t is the resonance production momentum transfer. We will present the details of 
the fitting procedure and the estimate of the theoretical uncertainties in the forthcoming 
paper [21]. Here, we give, for indication, the result with the following conditions: i) 
we adopt a set of harmonic oscillator wave functions with a common h.o. radius R = 
1/(0.4 GeV) = 2.5 GeV""*^ {"^{r) oc ea;p(— r^/2_R^)), which, we would say, is a widespread 
and empirically satisfactory recipe ii) we use 7 ~ 4 iii) we consider only low t data for the 
-ft'i(1270) decay, iv) we assume all the particles are on shell, except for the -ft'i(1270) — )■ Kp, 
v) we use the damping factor to he /3' = 3 GeV~^ (see the discussion on the damping factor 
in Appendix B). As a result, we find that the mixing angle 50° — 60° is well compatible 
with the data. Then, the partial wave amplitudes can be given, for example for 6k = 60°, 
as 

^l:i(1270)^K*7r ~ 1-6, ^Ki{1270)^K*7T ~ -0-2 

^!R'l(1400)^A'*7r ~ 3.1, y4^^(;^4Qo)_>.^,^ ~ 0.2 ^^^^ 

^I:i(i4oo)-^i<'p ~ —0.5, AKi{iioo)^Kp ~ 

4.3 Discussions on the phases and the scalar resonances 

4.3.1 The relative phases between different A'l — )• VP couplings 

The issue of the phases of the resonance couplings is very important in the present ap- 
proach, since the dependence on relies entirely on the phase of J'. And, indeed, one 
finds that changing the relative sign of the decays of i^i(1270) to Kp over K*^ channel 
would entirely change the prediction for A^. Then, we formulate the following observa- 
tions: 



• Phases between all the various resonances and decay channels into Knir are in 
principle measurable, and indeed are measured by ACCMOR collaboration. They 
can be also determined from our ^Pq model. From a theoretical point of view, one 
must not forget that, to obtain the full coupling sign of a quasi two-body channel, 
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one must take the product of the amphtude for the decay into the quasi two-body 
channel with the one of the isobar decay (e.g. K* — )■ Kn), in order to have the sign 
into the common final state Kirir, the only to make sense in comparing channels. 

However, the question remains far from trivial, because : 

• There were some misunderstandings in interpreting the ACCMOR data (e.g. the 
large D/S phase read by Gronau et al. does not in fact correspond to the D/S 
relative phase for the couplings of the 1400 ; the strong bump around 1400 in the 
D wave phase diagram. Fig. 13 in ref. [16], does not correspond to the D wave of 
the A'i(1400), which is very small ; it is a reflection of what happens in the S wave, 
since D wave phase is defined by reference to the S wave.). 

• Our model predicts real phases for all the couplings, which is also almost the case 
for the true i^— matrix predictions of ACCMOR collaboration. On the other hand, 
the data of ACCMOR show something different; so-called off-set phases, imaginary, 
and not predicted by the true i^'— matrix, are to be added to describe the data . The 
origin of these additional phases is unknown. It is very important to realise that 
that the solid lines of the histograms in [16] do not represent the true /T— matrix 
predictions, but include the ad hoc offset-phases. 

• One can test the soundness of our model by checking whether the predicted relative 
signs of couplings to K*7r and Kp agree with those shown by Daum et al.. In our 
study, we tried to establish the connection between the conventions of our model 
and those of ACCMOR collaboration. For the /A^ ratios in the common channel 
K*7r, the relation is trivial, and we find that there is agreement : A^ /A^ < for 
the J'S:i(1270), A^/A'^ > for the Ki{UOO) §. On the other hand, for the relative 
sign between the amplitudes of iri(1270) {K*n)s and iri(1270) {Kp)s, the 
conventions used by ACCMOR are not obvious, while this sign is crucial. In our 
study, we use the signs of ^Pq? but we also test different combinations of these 
relative signs, by allowing for additional phases 5. 

4.3.2 The controversial Ki(1270) — )■ Kutx decay through a scalar meson: Ki — )■ 

scalar + tt 

The PDG assigns a large branching ratio to this decay channel: i3r(i^i(1270) — )■ -ftrQ(1430)7r) = 
(28±4)%. It is extracted as all the branching ratios, from the ACCMOR data and analysis. 
However, this interpretation has been questioned. The original ACCMOR measurement 
shows indeed a clear, strongly coupled peak in the (scalar -|- vr) channel around the mass 
1270. However, it is not at all claimed that the scalar is i^'Q(1430) ; it is treated as a lower 
and much broader scalar meson (F ~ 600 MeV ; or could be a continuum (-ft''/r)s--„ave 

§In the latter case, the ACCMOR sign is deduced from the reanalysis by Babar [IT] 
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according to |25]. Indeed, our model predicts the decay to the J<'g(1430)7r channel to be 

of the order of 1%. Recently, the Belle collaboration has made a new branching measure- 
ment using the B — )■ J/iIj{iIj')Ki decay followed by Ki — )■ Kirir. What is most striking is 
that indeed. Belle finds Br{Ki{1270) i\:^(1430)7r) ~ 2% [E], as we predict, while not 
finding any new component in the Ki decay: the Br missing with respect to ACCMOR 
seems to be filled by an enlargement of Kp. Therefore, in our analysis, we do not include 
the -R'i(1270) — )■ ii'Q(1430)7r channel. We do not include either any other possible scalar 
in the presented results. However, to take into account the contradictory conclusions of 
ACCMOR we have kept in mind the possibility that there is some significant portion of 
the branching ratio carried by a very wide scalar meson, different from the i^Q(1430), 
such as the low lying state fi:(800). 

5 Determination of in the DDLR method 

In this section, we demonstrate how to determine the polarization parameter A-y from 
the actual experimental data using the maximum likelihood method. In particular, we 
introduce the DDLR method which was first applied in the r polarization measurement 
at the ALEPH experiment [S]. In the maximum likelihood method, knowing the 
dependence on the probability density function (PDF), the A^ closest to its true value can 
be obtained where the likelihood function (or equivalently, log-likelihood) given by the 
N sample of data takes its maximum value. In our case, the PDF, W, can be given as 
the decay width integrand normalized to unity (after due multiplication by the modulus 
squared of the Breit- Wigner) . Let us reiterate our statement that when one remains within 
the bump of the Ki resonance, the decay amplitude weakly depends on s = p^{Ki), and 
one can set s = ml-^ in their expression, i.e. in the jT's, which we assume therefrom. 
Thus, using Eq. (ITT]) , we find 



W^(si3, S23, cos 9) = /(si3, S23, COS 9) + \^g{si3, S23, cos 9) 



(32) 



where 



/(S13,S23,C0S6') 



1 



\J\\l + cosH) 



g{si3,S23,COs9) 



4/ 

1 



Im[n- {Jx J*)] cos 9 



(33) 
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where /, g are normalised relatively to the measure dsi3ds23dcos9. 

Then, the likelihood function for the TV events of data can be given as 



N 




(34) 



i=l 



14 



where i indicates the kinematic variable of each event. The true value of should 
maximize this function, namely should be the solution of the following equation: 

dC 



(35) 



Now, we explain the DDLR method [5]. The next procedure to look for the value of 
A^ in our problem is usually to use the known distribution of / and g functions and fit 
the value of A^ so as to maximize the likelihood function. It should be noted that this is 
not a very simple task, especially since the / and g are complicated functions as shown 
in Eq. (15^ and the Appendix A. In [3], it is pointed out that when the PDF depends 
on the parameter, which we are interested in, only linearly, one can reduce such a multi- 
dimensional fit to an one-dimensional one using a single variable u which is defined as 
follows: 

/ /)N ^(S13,S23,C0S^) 

W(-Sl3, S23, COS 9) = — — . (36) 

/(Sl3,S23,COS^) 

This can be proved simply by writing down the log-likelihood of our problem: 

N 

ln£ = ^ln [1 + A^ Cl'(s^3, S23, cos ^*)] + terms independent of A^. (37) 

i=l 

where A^ does not depend on / and g separately but only on their ratio u. This demon- 
strates that only the oj distribution is needed to extract A^. 

Another thing that is pointed out in [5] is that the polarization parameter is often 
determined only by using the angular distribution, however, the sensitivity to it can be 
further improved by considering all the kinematic information, such as the Dalitz variable 
distribution. Therefore, we use the dependence of A^ not only on cos6' but also on S13 and 
S23 in this work. Considering the fact that /, g and u have very complicated dependences 
on these kinematic variables, the reduction to the one-dimensional fit achieved by using 
the variable u is very efficient for the data analysis as shown in the following. 

It must be underlined that, in the present case, in contrast with the initial DDLR 
problem, r — )■ ttz/, u is not a purely kinematic variable, it depends itself on the theoretical 
model, as it is the case for t ^ aiu. Then the method to obtain the distribution in u is as 
follows. First, following the standard MC method, we generate the faked events according 
to the PDF. Then, we compute the u value for each event. In this way, we obtain the 
omega distribution N^'-^{uj) according to the PDF. 

We show examples of the u distribution generated by the Monte Carlo (MC) simulation 
in Fig. El 

Now we explain how to extract the value of A^ as well as its statistical error from a 
given CO distribution. We will we present our sensitivity study result in section O Since 
the use of the u variable reduces our fit to a one-dimensional one, A^ is obtained simply 
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Figure 2: The simulated w-distribution for = +1 (red) and A^ = —1 (blue). The 
polarization parameter A^ can be determined from the difference between these two dis- 
tributions (see the footnote for more details). 



as a solution to the following equation: 

oX^ 1 + X^u 

Of course, one could solve the equation fl38l) by successive searches. However, we can 
provide explicit expressions for A^. One sees easily that the normalised distribution in w, 
W'{(jj), can be written as : 

W\u) = (t){uj){l + X^u) (39) 

where (j){uj) is an even function of w, since as can be seen from Eq. (155]) /(sis, S23, cosO) 
is an even function of cos^ while Ci;(si3, -5235 cos6') is an odd one. 

Then, one can easily demonstrate by integration over the interval —1 < uo < 1 that A^ 
can be expressed as ratios of odd over even momenta: 

Therefore, the expression obtained by DDLR for small A-y seems exact. 

Similarly to Eq. f l38p . one can also obtain the statistical error to the given value of A^ 

as: 

< = rr- (41) 
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Thus, once the u distribution is obtained experimentally, the Eqs. (138|) or ( 140|) and (l4T]l 
immediately provide the values of and . 

6 Future prospects for the polarization measurement 

In this section, we discuss the sensitivity of the future experiments, namely the SuperB 
factories and the LHCb to A^, using the B — )■ i^i (1270)7 ~^ (A' 7171)7. We also discuss 
the advantages and disadvantages of our method compared to the other methods of the 
polarization measurement using the other processes, such as -B — > K*e^e~ , Bd K*'y 
and 5s — )■ 07- 

6.1 The sensitivity study of the polarization measurement with 
B 7^1(1270)7 in the DDLR method 

In this section, we perform a Monte Carlo simulation in order to estimate the sensitivity 
of the future experiments to the polarization parameter A^ using the DDLR method. 
Following the procedure described in Section 5, we first generate the events (10'^ and 10^ 
events as examples) for a given value of A^ and then estimate the expected statistical error, 
(Ja^. Here, we use the "ideal" Monte Carlo simulation, i.e. detector and background effects 
are not taken into account. In order to generate the events as well as to compute the ou 
distribution, we use the input hadronic parameters as given in the Section 4, taking into 
account the form factor momentum transfer dependence (discussed in Appendix B). These 
parameters include the experimentally measured isobar widths, the ^Pq model parameters 
(the meson wave function radii, the quark-pair-creation constant, damping factor) and 
the phenomenological Ki mixing angle. In the Table [H we present our result in the case 
of the SM, i.e. A^ = 1. One can see from the table, that for 10^ events the error on A^ 
is smaller than 0.1. We found that the errors do not change much for different values 
of A^. We found that the u distributions for the K~^7T~^7r~ and K^7t^it~, and K^ir^ir^ 
and K~^7r~n^ are the same. Then, it should be pointed out an advantage of using the 
a;-variable: all the channels corresponding to the same PDF can be merged altogether. 

^In the real data, one must consider the systematic errors coming from detector effect etc and perform 
a fit instead of using these simple formulae. There is one subtlety for that case. For each event, 
the photon should have the polarization either left- or right-handed. Thus, in MC, we produce the uj 
distribution with purely left- and right-handed PDF. Then, the total uj distribution of the experimental 
data is expected to be a linear combination of these two distribution with a ratio of e: 

N^^'P-iiu) = eN^^iuj) + (1 - e)Nj^^{io) 

with e = — The N is the number of event in the experimental measurement. We show an example of 
the UJ distribution of = —1 (red) and = +1 (blue) in Fig. [2] As seen in this equation, the can 
be determined from the difference between these two distributions. 
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That means that one can compute w- variable for each event and build a single histogram, 
which can increase the statistical significance. 

In the above, we use the full decay distribution, not only on the information of the 
angular part but also the information of the invariant mass of the hadronic system. In 
the original DDLR paper [5], it was pointed out, using an average decay distribution in 
place of a full decay distribution for each set of invariant masses results in a decrease of 
the sensitivity. In order to test this, we also produce the u distribution including only the 
cos6' dependence, i.e. integrated over the Dalitz plot, and compute cta^. We found that 
the inclusion of the full Dalitz information can indeed improve the sensitivity by typically 
a factor of two comparing to the angular fit. 

Up to now, we have not considered the systematic errors coming from the hadronic 
parameters. We must reiterate that our hadronic model applied in the above analysis 
is approximate; it depends on basic assumptions like the non relativistic approximations 
inherent to the quark models. It depends also on parameters, some of them being internal 
to the full quark model, like the meson radii, and one being purely phenomenological, the 
mixing angle 6ki (we must note that there exists a correlation between the mixing angle, 
extracted from the data, and the chosen set of meson radii). It depends also on the set of 
experimental data which we claim to describe by such models as discussed in Section 4. It 
is then a difficult question to evaluate the uncertainties of our results ; we do not claim to 
discuss this point precisely in this paper but we intend to do this in another publication. 

Finally, we would like to give a rough estimate for the event numbers expected by the 
future experiments, namely the SuperB factories and the LHCb. Taking the exclusive 
branching fraction Br{B~^ — K/' (1270)7) = 4.3 x 10~^ and assuming that the decays 
Ki — > Kim are by K*7r (16%) and Kp (42%) channels, we obtain the observable branching 
fraction of fir(fi+ {K+7r-7i+)K,(i270)l) = 4.3xl0-^x(0.16*4/9+0.42*l/6) ~ 0.6x10"^ 
and Br{B+ ^ (i^°vr+7r°)K^(i27o)7) = 4.3 x 10"^ x (2 * 0.16 * 2/9 + 0.42 * 1/3) x 1/3 ~ 
0.3 X 10~^ (here the last factor 1/3 comes from the fact that K'^ is observed as tt+tt" 
from the Ks decay). In order to get a more realistic estimation of the required number of 
signal events in the future experiments, we take the total efficiency of the reconstruction 
and selection to be of the order of 0.1% as in the case of i? — > K*'~f and Bg — )■ 07 at 
the LHCb experiment [2S] and of the order of 1% at B factories Then, we obtain the 
yield of the nominal data taking to be of the order of 5 x 10"^ of B'^ {K~^'n'^'^~^)Ki{i27o)'l 
and 2.5 x 10^ of B~^ — )■ {K^h~^7i'^)ki(i27o)1 signal events for the accumulated luminosity 
of 2 fb~^ at LHCb. The estimated annual yield at SuperB factories with 2 ab~^ of 
integrated luminosity is of the order of 1 x 10^ and 0.5 x 10'^ of signal events of B~^ — >• 
{K~^71~'7i~^)ki{i27o)'1 and B'^ — )■ (-ft'°vr+7r°)/^^(i27o)7, respectively. Thus, the event sample, 
10'^ and 10^, studied in Table 1, roughly corresponds to the annual expected events of 
SuperB and LHCb, respectively. It should be noted that the decay modes including a 
neutral particle is difficult to be studied in LHCb, i.e. LHCb may well study the first 
decay channel in the Table 1 while SuperB can study all of them reasonably well. 
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/ J- J- • J- • 1 \ 

(statistical error) 


events 


1\T T m4 

events 


B+ (i^+7r-7r+)i^,(i270)7 


± 0.18 


± 0.06 


B+ ^ (X°7r+7r0)^,(i270)7 


± 0.12 


± 0.04 


50 ^ (/i07r+7r-);,,(i270)7 


± 0.18 


± 0.06 


50 ^ (i^+vr-7rO)A',(i270)7 


± 0.12 


± 0.04 



Table 1: Sensitivity study of the polarization measurement with B — )■ 1270)7 iii 
the DDLR method. Our estimates of the statistical errors to in the case of SM (i.e. 

= +1) is shown in this table. The event sample, 10^ and lO'^, roughly corresponds to the 
annual expected events of SuperB and LHCb, respectively. The hadronic parameters used 
to obtain this result are given in section 4. The systematic error due to the uncertainties 
from these hadronic parameters is not included and has to be carefully studied. 



6.2 Comparison to the other methods 

In this subsection we compare the precision of the photon polarization measurement, using 
various direct and indirect methods. 



6.2.1 Comparison with the up-down asymmetry of GGPR 

One of the direct methods of the photon polarization determination methods, proposed 
by Gronau et al. [3], is to study the angular distribution in the B — t- PiP2-f37 decay 
and extract the polarization parameter A^ from the angular correlations among the final 
hadronic decay products P1P2P3. An observable called "up-down" asymmetry is intro- 
duced: 

L\^cos^^ I dsds^,ds,s\J\^ ^ ^ 



representing the asymmetry between the measured number of signal events with the pho- 
tons emitted above and below the P1P2P3 decay plane in the Ki reference frame. Having 
the theoretical prediction of JT, one can determine A^. 

Our conclusion, identical to the one for the angular fit, is that the statistical error on 



A-y is about twice the one in our method. 



6.2.2 Comparison with B K*t 



From the analysis of the angular distributions of the four-body final state in the B^ — )■ 
i^*'^(— >■ K~TT^)l'^ l~ decay in the low invariant mass region one can study various 
observables that involve different combinations of K* spin amplitudes 



Working in the transversity basis A^_l = ^^^j^^ and = ^^^^^ 1 one the most 
promising observable, that has a small impact from the theoretical uncertainties, is the 
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transverse asymmetry defined as 



(2) ^ \M±\^-\M\i\^ _ MrMI + M^Ml 



\M±\^ + \Mii\^ \Mr\^ + \Ml\^ 



SM 



Note that we assume that M.l/r at the low I'^l invariant mass regions can be identified 
with the decay amphtudes of 6 — s^iljr and are related to our polarization parameter as 

~ \Mr\^-\Ml? 



|2 



The new analysis of the B — > K*e^e~ decay mode by the LHCb collaboration 
shows that one can expect an annual signal yield of 200 to 250 events for 2 fb~^ in this 
energy region. With this number, it is found that the LHCb can reach a precision of 
a{Aj, ) around 0.2 corresponding to the statistical error on (j{M.l/ M.^) to be of the 
order of 0.1 [M]. 

It should be noticed that this method allows the direct measurement of the ratio 
X = \M.l/ M.r\, while our polarization parameter is sensitive only to the amplitude 
ratio square, x^. Therefore, the errors of these two methods are to be compared using the 
following equation: 

which shows that the sensitivity depends on the value of x. We should immediately 
notice that for verifying the SM value, x ^ 0, the method accessible to x is much more 
advantageous than the one to x^: our is in fact insensitive to the SM point (requiring 
an infinitesimal error). We plot Eq. HH in Fig. [31 Let us look at the horizontal line of 
ax = 0.1, expected error on x with the B — K*e~^e~ measurement. One can see that our 
method becomes more advantageous above x ~ 0.3, where the same sensitivity to x can 
be achieved with a larger error to A^ i.e. cr^^ > 0.1. 

6.2.3 Comparison with the methods invoking CP asymmetries 

An indirect method to measure the photon polarization is to study the time-dependent 
CP asymmetry in the neutral Bq {q = d, s) mesons. For the generic radiative decay 
of the neutral i?q-meson into any hadronic self-conjug ate state M^^, Bq{t) M^^-f, 
neglecting direct CP violation and the small width difference between the two S-mesons, 
a CP- asymmetry is given by 



Acp{t) = ^ sin(2?/') sin((^M - - <Pr) sin(Amt) (45) 

where ^(= ±1) is the CP eigenvalue of M*"^, sm.{2ip) = parametrizes the 

relative amount of left- and right-polarized photons, 0l,_r = sin^^ ( ^[a4l ) 
ative CP-odd weak phases in the b s'-f process and (pM is one in the Bg — Bg mixing. 
These phases are 4>l/r = 0, 4>d = 2/3, </)s — in SM. The smallness of the right-handed 
amplitude in SM, J^l/J^r ^ 2ms/mb, predicts Acp(t) ~ 0. We should emphasise that 
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Figure 3: Comparison of the sensitivity of the two methods: the one directly determining 
X = \M.l/ M.ji\ and the other one determining such as our (see Eq. (144|) ). One can 
see that when we assume the same errors for the both methods, a better significance can 
be obtained with the later method for x > 0.3. 

Acp(t) measures the combination of x = \A4l/-Mr\ and the CP violating phases 4>m,l,r 
but not separately. Thus, the value of x can be obtained from this measurement, only by 
having the value of the CP violating phases in the b ^ S'j as well as the Bq mixing. 

Current world average for the asymmetry in the — )■ Ks'ir^'y process is Scp{B — >■ 
Ks'^'^'y) = —0.15 ±0.20 [T], which is expected to be improved by the SuperB factory; error 
down to 2% [20]. The LHCb experiment is going to measure Bg — )■ 07 process. Based 
on the MC simulation for 2 fb^^, it is claimed in [26] that LHCb will be able to measure 
X with the accuracy of — 0.1. Therefore, similar to the case of B ^ K*e^e~, our 
method using can be more sensitive to x above x ~ 0.3 (see Fig. |3]). Again, it should 
be emphasised that although an observation Acp{t) 7^ in this method immediately 
indicates the existence of new physics, a quantitative determination of x is not possible 
unless we fix the new physics model, namely the CP violating phases in 6 — > S7 as well 
as Bq mixing. 

7 Conclusions 

We investigated the method to determine the photon polarization of the 6 — j- 57 process 
using the decay channel B — )■ Ki^y KinT'-f, which was originally proposed by Gronau 
et al [3]. In this paper, we propose a new variable, u, to determine the polarization 
parameter A^. This variable was firstly applied in the r polarization measurement in the 
ALEPH experiment [5]. The use of u significantly simplifies the experimental analysis 
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and as a result, it allows us to include not only the angular dependence of the polarization 
parameter, represented as the up-down asymmetry in [3], but also the three body Dalitz 
variable dependence to the fit. We found that when the data is analysed by using w, 
the statistical error in the polarization parameter can be reduced by a factor of 2, 
comparing to the case of the up-down asymmetry. 

In order to evaluate the systematic error, a sufficiently accurate modelling of the 
hadronic decays of Ki — )■ A'tttt is required. Having the recent observation of the Belle 
collaboration |1] implying the dominance of the B — )■ i^i (1270)7 channel over B — )■ 
i^i(1400)7, we investigated the hadronic decay of ii'i(1270) — > Kinr in some detail. 
We first have derived the basic hadronic parameters required in our analysis. These 
parameters can, in principle, be determined by the experimental measurements of the 
J'ri(1270) —7- Knir decay. On the other hand, although the outstanding ACCMOR ex- 
periment had provided an extensive study of this decay, we found that the information 
one can extract from it is not accurate enough. We described some of the problems en- 
countered in our analysis, which include the strong phase between different intermediate 
resonance states and the controversial -R'i(1270) — )■ Kmr through scalar mesons. Interest- 
ingly, the latter problem has been studied by the Belle collaboration recently [IB] , with a 
small result contradicting to the PDG number and in an agreement with our prediction 
by the ^Pq model. Their results will provide a great help for our future study of the Ki 
decay modes. 

Being unable to obtain the hadronic parameters from the fundamental theory, we 
resorted to combine experimental data and phenomenological models. Practically, com- 
bining the experimental results of the partial wave analysis of the Ki decays and the 
predictions of the ^Pq quark-pair-creation model, we found that the Ki mixing angle 
between 50° to 60° is well compatible with the experimental data. Nevertheless, an eval- 
uation of the theoretical uncertainties requires much more detailed discussions, which we 
will present in the forthcoming paper. 

Appendix A: The Ki — ^ Ktttt form factors 

In this appendix, we derive the quasi-two-body decay amplitude J^Jp.p.^p^ given in Sec- 
tion 13.21 For the computation of this amplitude, we take into account the vector meson 
resonance width effect assuming the Breit-Wigner form, thus 

Mjp^P^^P^ ^ M{Ki ^ VPk)M{V ^ P,P,)BWv{sij) (46) 

The decay amplitude of the axial-vector Ki to a vector {V) and a pseudoscalar (P^) meson 
can be expressed in the following Lorentz invariant form: 

{V{pv,e^''^)PkiPk))\AHA\K,{pK,,e^'^^^)) = ef^^T^-e^y> , T^- = fyg^^ ^hyv'yV\, (47) 
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where fv and gv are the form factors. The amphtude of the subsequent decay V to two 
pseudoscalar mesons Pi and Pj can be parametrized in terms of one vector-pseudoscalar 
couphng gvPiPj- 



{P,{Pi)P,{p,)\AHy\V{pv,e^''^)) = gyp^p^e(y\p^-p.Y 



(48) 



Using these form factors, we can obtain 



MY, 



{Pipj)pk 



{Pi ■ ^Kja^j + {pj ■ eK,)k 



(49) 



where 




gvp,p^BWv{sij)[fv + hvVs{Ei - Ej) - Aij] 
gvp.p,^'^v{sij)[-fv + hv^fs{Ei - Ej) - A^,] 



(50) 



with A,,- = 




[fv + hvyfs{E, + Ej)]. Note E,, = (s - s^u + m^)/{2y^). 



Finally, using these afj and bfj functions, we obtain the Ki — P1P2-P3 amplitude as: 



Appendix B: Damping factor 

In this appendix, we discuss the necessity of introducing the cutoff in the coupling vertices. 
When we compute the transition rates, we must take into account, in principle, the 
widths of the initial or final resonances; this is especially crucial for the transition rate 
of i^'i(1270) — Kp, which is large, although it would be kinematically almost forbidden 
at the nominal values of the masses. A well-known and simple way to take widths into 
account is by integrating over the off-shell "masses", p^, with the weight of the Breit- 
Wigner's. However, it is then found that the integrals will diverge for P or D waves, due 
to the /c^' factors, where k is the decay momentum, if the coefficients are taken to remain 
constant. Of course, the reactions will in general provide natural limits of integration, 
for instance the spectrum studied by ACCMOR stops at Mxmr = 1.6 GeV, but even 
that cut would give exceedingly large P ot D wave contributions. In fact it seems that 
various indications hint to the necessity of a strong dynamical cutoff, or "damping factor" , 
affecting for instance the Breit-Wigner shape (e.g. accurate studies of A (1236) [3T] or 
i^*(890), see ref. [32]); the prototype of which are the Blatt-Weisskopf factors. The need 



■M{Kl — !■ P1P2P3) — Cl32A^(PiP3)P2 + C23lA^(P2P3)Pi + Cl23A^(PiP2)P3 

= {Pi ■ ^kJCi - {p2 ■ eKjC2 



(51) 



where 



Cl — Ci32(ai3 — &13) — C231&23 + Cl23'2l2 
C2 = Ci32&]^ — C23i(a^ — 6^) — Ci23&]^ 



(52) 
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for it is also shown by calculations of hadronic loops in the ^Pq model [33]. One obtains 
a natural damping factor through the gaussian factors e~^^'^: 

oc (3-aA;2)g-/3fc2^ oc afcV^'^' (53) 

but /3 ~ 0.3 GeV~^ is found there to be quite too small. Following ref. |33], we introduce 
the empirical Gaussian cutoff exp(— /S'/c^), renormalizing (3 ^ (3 + (3' and we fix (3' ?a 
3 GeV-2. 

Now, to study the Ki decays, we integrate the squared amplitudes over the Ki and 
vector meson resonance states {K* /p) invariant masses within the whole kinematic allowed 
region in ACCMOR, [1.0; 1.6] GeV; the integration on Ki invariant mass does not depend 
too much on these limits for S waves. It is not the case for the D waves, but once the 
damping factor is introduced, the D/S ratio becomes stable. The isobar {K*/p) decay 
does not depend much on the damping factor. 
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